PASJ: Publ. Astron. Soc. Japan , 1-??, 
(c) 2013. Astronomical Society of Japan. 



X-Ray Spectrum of a Peculiar Supernova Remnant G 359.1—0.5 
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Abstract 

We present the Suzaku results of a supernova remnant (SNR), G 359.1—0.5 in the direction of the 
Galactic center region. From the SNR, we find prominent K-shell lines of highly ionized Si and S ions, 
together with unusual structures at 2.5-3.0 and 3.1-3.6 keV. No canonical SNR plasma model, in either 
ionization equilibrium or under-ionization, can explain the structures. The energies and shapes of the 
structures are similar to those of the radiative transitions of free electrons to the K-shell of He-like Si and 
S ions (radiative recombination continuum: RRC). The presence of the strong RRC structures indicates 
that the plasma is in over- ionization. In fact, the observed spectrum is well fitted with an over- ionized 
plasma model. The best-fit electron temperature of 0.29 keV is far smaller than the ionization temperature 
of 0.77 keV, which means that G 359.1—0.5 is in extreme condition of over-ionization. We report some 
cautions on the physical parameters, and comment possible origins for the over- ionized plasma. 

Key words: Galaxy: center — ISM: abundances — ISM: individual (G 359. 1—0.5) — ISM: supernova 
remnants — X-rays: ISM 



1. Introduction 

G 359.1—0.5 was discovered by Downes et al. (1979) as 
a shell-like object with the 4.875 GHz radio band, and 
was suggested to be a supernova remnant (SNR). Reich & 
Fuerst (1984) observed G 359.1-0.5 with the 2.695 GHz 
and 4.75 GHz radio bands, and found a complete shell 
structure. Together with the spectral index (a = —0.37) 
and the liner polarization, they confirmed this object to 
be an SNR. Uchida et al. (1992a) found a ^^CO molecular 
ring that is concentric with the shell of G 359.1—0.5. The 
molecular ring has high radial velocities between —60 and 
— 190 km s~^. Furthermore, Uchida, Morris, and Yusef- 
Zadeh (1992b) found a ring structure of G 359.1—0.5 in 
the 1.42 GHz band, surrounded with the Hi emission ring 
having large negative velocities of —75 to —190 km s~^. 
These large velocity differences of about 100 km s~^ would 
be due to the shear motions in the Galactic center (GC) 
region, and hence G 359.1—0.5 is likely located in the GC 
region. 

The non-thermal radio filament, named the "Snake" 
(e.g. Gray et al. 1991; Gray et al. 1995), is located at 
the projected position of G 359.1—0.5. Along the shell of 
G 359.1—0.5 and perpendicular to the Snake, Lazendic et 
al. (2002) discovered shocked molecular hydrogen (H2). 
The intriguing fact is that an OH (1720 MHz) maser 
source, a signature of an SNR/molecular cloud interac- 
tion, is located in the H2 cloud at the projected inter- 
section of G 359.1—0.5 with the Snake (Yusef-Zadeh et al. 
1995). Thus G 359.1—0.5 and the Snake may be physically 
interacting, although it could be an artifact of superposi- 



tion of the disparate components along the line of sight. 
These radio morphologies place G 359.1—0.5 as one of the 
enigmatic SNR possibly related to the extreme environ- 
ment of the GC. 

The X-ray emission inside the radio shell was first dis- 
covered with ROSAT (Egger & Sun 1998). The ASCA ob- 
servation also found the center-filled X-rays and claimed 
prominent K-shell emission lines from He-like ions of Si 
and those from H-like ions (Lya) of S (Bamba et al. 2000). 
From the radio and X-ray morphologies, G 359.1—0.5 is 
likely to be a mixed- morphology SNR (Rho & Petre 1998), 
although the X-ray image data were poor for this concrete 
statement. 

The line features of ASCA were very peculiar, because 
S should be ionized far higher than Si, although S and 
Si are synthesized in the same layer of a progenitor star. 
In fact, the X-ray spectrum required a two-component 
collisional ionization equilibrium (CIE) plasma of 0.6 and 
4.4 keV temperatures and an extreme over- abundance of S 
(> 38 solar) in the higher temperature plasma. Thus the 
X-ray spectrum of G 359.1—0.5 is also enigmatic. These 
peculiar features in the radio morphologies and in the X- 
ray spectrum may be related with each other and may 
require a new scenario for the structure and evolution of 
the SNR. 

We observed G 359.1—0.5 with Suzaku (Mitsuda et al. 
2007). With the excellent energy resolution, the large 
effective area, and the low background for diffuse emis- 
sion, the X-ray Imaging Spectrometer (XIS; Koyama et 
al. 2007b) provide the first opportunity to entangle the 
peculiarity of G 359.1—0.5. In this paper, we report the 
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high quality X-ray spectra and detailed analyses, and dis- 
cuss the nature of G 359.1—0.5. All of the errors are at 
the 1(7 confidence level throughout this paper, unless oth- 
erwise mentioned. 

2. Observations and Data Reduction 

The G 359.1—0.5 region was observed as a part of 
the Galactic Center Survey Project of Suzaku on 2008 
September 14-16. We also used its vicinity to obtain a 
background spectrum and to examine surroundings of the 
G 359.1-0.5 (Bamba et al. 2009). The observation logs 
are given in table 1. 

We used the XIS data for the imaging and spectral 
analyses. Details on the XIS are found in Koyama et 
al. (2007b). Further information of XIS on the calibra- 
tion, qualification, and so on, is found in separate papers 
(e.g., Uchiyama et al. 2009). Although the XIS consists 
of four CCD cameras with three front-illuminated (FI) 
CCDs (XIS 0, XIS 2, and XIS 3) and one back-illuminated 
(BI) CCD (XISl), we did not use one of the FI CCDs 
(XIS 2) because of the loss of the proper function since 
2006 November. 

The data processing and screening were made using the 
processing version 2.2.11.22, xispi software version 2008- 
06-06, and the calibration database version 2009-07-17. 
After the processing and screening under the standard 
criteria^, the exposure time were 57.7, 70.5, and 72.5 ks 
for the G 359.1—0.5 and two vicinity fields, respectively. 

3. Analysis 

We use HEAsoft version 6.6.1 and SPEX (Kaastra et 
al. 1996) version 2.02.02 in the imaging and spectral anal- 
yses. The Non- X-ray Background (NXB) is made using 
xisnxbgen (Tawa et al. 2008) and subtracted from the 
raw data. XIS response and auxiliary files are created us- 
ing xisrmfgen and xissimarf gen (Ishisaki et al. 2007), 
respectively. Since we use SPEX for the spectral fitting, 
we make format conversion using traf o version 1.02. 

3.1. Image 

Referring to the X-ray spectrum of G 359.1—0.5 by 
Bamba et al. (2000), we show the X-ray images in the 
(a) 1.5-2.3 keV and (b) 2.5-3.0 keV bands. After the 
NXB subtraction and the subsequent corrections of the 
exposures and the vignetting effects, all the XIS images 
are merged to increase statistics. We shift the Suzaku co- 
ordinates referring to the XMM-Newton sources (Watson 
et al. 2009), because the Suzaku nominal position error 
is ^ 20^' (Uchiyama et al. 2008), on the other hand, that 
of the XMM-Newton catalog is - 2". The name of the 
reference stars detected with XMM-Newton and Suzkau, 
their positions (/, b) and differences (A/, Ab) are listed in 
table 2. After these fine-shifts of the Suzaku coordinates, 
we make the images in figures la and lb. The figures 
are binned with 8f'3 x 8f'3 and smoothed with a Gaussian 

^ http: / /heasarc. nasa.gov/docs / suzaku/processing/criteria_xis.hti 
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Fig. 1. XIS images in the (a) 1.5-2.3 keV and (b) 2.5-3.0 keV 
bands. The source and background regions are shown as a 
sohd eUipse and a dotted pentagonal (white), respectively. 
The radio map by the VLA 90 cm observation is superimposed 
in green contours. 

kernel of <j = 0.^70. The radio map by the VLA 90 cm ob- 
servation (LaRosa et al. 2000) is superimposed with green 
contours. 

In figures la and lb, we find diffuse X-rays inside the 
radio shell of G 359.1—0.5 (green contours). No excess 
X-ray relative to the surrounding region is found on the 
radio shell. This is the first convincing evidence that 
G 359.1—0.5 is a mixed-morphology SNR. We designate 
the "source region" by the solid ellipse centered on 
6) = (+359.° 10, — 0.°49) with the semimajor and semimi- 
nor axes of 7'9 and 4f9, respectively. The bright emis- 
sion in the lower corner of figure la is a part of another 
SNR G 359.0-0.9 (Bamba et al. 2000; Bamba et al. 2009). 
Since the present observations cover only a small portion 
of the SNR, we do not treat this SNR in this paper. 

3.2. Spectra 

The X-ray spectra of G 359.1—0.5 are made from the 
source region. Since this SNR is located in the direction 
of the GC, the Galactic Center Diffuse X-ray (GCDX; 
e.g., Koyama et al. 2007a) is a major background for the 
spectrum. We therefore select the dotted pentagonal re- 
gion as shown in figures la and lb, where no X-ray excess 
is found (hereafter, the background region). The mean 
Galactic latitude is almost the same as that of the source 
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Table 1. 


Observation loj 
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Observation ID 


— c J 

Start date 


R.A. 


Dec. 


i — 

rL/xposiire 






(J2000.0) 


(ks) 


503012010 


2008-09-14 


266.°2992 


-29.°9408 


57.7 


502016010 


2008-03-02 


266.°2249 


-30.°1096 


70.5 


502017010 


2008-03-06 


266.°4705 


-30.°0867 


72.5 


* Effective exposure 


of the screened 


XIS data. 







Table 2. Reference sources used in the position adjustments of the Suzaku data. 



Observation ID 


Reference 


source 


Source position 


Source position 


Difference 


Suzaku 


XMM-Newton 


XMM-Newton 


Suzaku 














b 




b 


Al 


Ab 


503012010 


J174539.4- 


-300139 


359.°0726 


-0.°5756 


359.°0733 


-0.°5765 


0.°0007 


-0.°0008 


502016010 


J174458.3- 


-300655 


358.°9205 


-0.°4951 


358.°9207 


-0.°4952 


0.°0002 


-0.°0001 


502017010 


J174558.0- 


-295738 


359.°1648 


-0.°5982 


359.°1674 


-0.°6003 


0.°0026 


-0.°0021 



region, which may compensate the strong flux dependence 
of GCDX on the Galactic latitude. The source and back- 
ground spectra are made by subtracting the NXB. After 
the correction of the differences of the effective areas and 
the effective exposures including the vignetting effects for 
both the source and background region, we subtract the 
background spectra from the source spectra. To improve 
the photon statistics, the background-subtracted source 
spectra of the two FIs (XIS and XIS 3) are merged be- 
cause the response functions are almost identical. The 
background-subtracted FI spectrum of G 359.1—0.5 (the 
source region) is shown in figure 2a. 

We see no significant X-rays in the hard X-ray band 
above 4 keV, and hence fit the FI and BI spectra in the 
1.0-4.0 keV band simultaneously, although we show only 
the FI spectrum in the figures for brevity. 
3.2.1. ColUsional Ionization Equilibrium plasma 

Using the code of SPEX, we fit the spectra with the 
model of an optically thin thermal plasma in CIE af- 
fected by photoelectric absorption with the cross-sections 
of Morrison & McCammon (1983) (1-CIE fit). The abun- 
dances of Si and S {Zsi and Zs) are free parameters, but 
those of the other elements are fixed at the solar values 
(Fodders et al. 2009). The best-fit 1-CIE model is given 
in figure 2a with the solid line, while the data residuals 
are shown in figure 2b. Residuals are found around the 
energies of 2.5-3.0 and 3.1-3.6 keV, although the latter is 
marginal. Therefore, this model is rejected with the large 
X^/ degrees of freedom (d.o.f.) of 585/191. 

We then fit the spectra with a two-component CIE 
model with the common absorption for both the CIE com- 
ponents (2-CIE fit). The abundances of Si and S are free 
parameters but are Hnked between the two components. 
The abundances of the other elements are fixed at the solar 
values. The data residuals from the best-fit 2-CIE model 
are shown in figure 2c. The residuals around 2.5-3.0 keV 
are still found even by this model with the /d.o.f. of 
365/189. We thus try the 2-CIE model fitting but with 
unlinked Si and S abundances. This is the same model ac- 
cepted for the ASCA spectrum (Bamba et al. 2000). We, 



however, find the unacceptable fitting with x^/d.o.f. of 
360/187, remaining the significant residuals as is shown 
in figure 2d. 

The best-fit values of the cooler and hotter plasma tem- 
peratures are 0.20 oi ^-^ ±0.1 keV, respectively, 
which are inconsistent with those reported by Bamba et 
al. (2000). This inconsistency would be due to the differ- 
ent sensitivity in the energy band below Mg-Lya energy 
(< 1.5 keV) between the Suzaku XIS and ASCA CIS (Gas 
Imaging Spectrometers; Ohashi et al. 1996). Although 
Bamba et al. (2000) fitted the spectra in the energy band 
of 1-1.5 keV, the sensitivity of this energy band is very 
small. The XIS, on the other hand, has higher sensitivity 
in the same energy band. We therefore fit the XIS spectra 
with the best-fit model and values reported by Bamba et 
al. (2000), and find large residuals at the low energy and 
around 2.5-3.0 keV with xV^.o.f. of 821/193. 

3.2.2. Under-ionized plasma 

Since almost all the plasma of young SNRs are in under- 
ionization, we fit the spectra of G 359.1—0.5 with one or 
two-component under-ionized plasma model affected by 
photoelectric absorption using the SPEX code. The abun- 
dances of Si and S are free parameters, while the abun- 
dances of the other elements are fixed at the solar val- 
ues. These model are rejected with the large x^/d.o.f. of 
585/190 and 365/188. In the both cases, the best-fit ion- 
ization temperatures {kTz) are the same as those of the 
electron temperature (kTe). These fittings are essentially 
the same as those of the CIE model. We therefore do not 
proceed further analysis on the under-ionized plasma. 

3.2.3. Over-ionized plasma 

The energies of the residuals found in the 1-CIE fit cor- 
respond to the K-shell binding potentials of the He-like 
Si (2439 eV) and He-like S (3225 eV). These structures 
are similar to the radiative transitions of free electrons to 
the K-shell of ions (radiative recombination continuum: 
RRC) found in IC443 and W49B (Yamaguchi et al. 2009; 
Ozawa et al. 2009). From the energies of the structures, 
these are thought to be the RRC for the transitions of free 
electrons to the K-shell of He-like Si and S ions. We there- 
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Fig. 2. (a) Background-subtracted FI spectrum of 
G 359.1-0.5 with the best-fit 1-CIE plasma model. (b) 
Residuals from the best-fit 1-CIE plasma model. (c) 
Residuals from the best-fit 2-CIE model with the common 
abundances, (d) Same as (c), but the Si and S abundances 
are unlinked between the two components. 

fore examine if the RRC are real or not, by adding the 
RRC of He-like ions of Si and S to a 1-CIE plasma model. 
For consistency, we also add the Lya lines of Si and S 
with the fixed energy of the theoretical values at 2006 and 
2622 eV, respectively (hereafter. Model A). The electron 
temperature values for the RRC of Si and S are common 
free parameters. Then the x^/d.o.f. value is dramati- 
cally improved to 224/186, where the best-fit spectrum of 
Model A are shown in figure 3. The best-fit CIE temper- 
ature and abundances are kTciE of 0.62 06 ^^^^ ^si of 
4.1 1q"5 solar and Zs of 4.2 solar, respectively, while 
the best-fit RRC temperature is kT^ of 0.31 to^^l keV. 

As is shown in figure 3, we find the RRC of He-like 
Si ions mainly contribute to the structure, of which the 
ASCA observation claimed the Lya lines of S ions (Bamba 
et al. 2000). Unlike IC443 and W49B, the fluxes of the 
RRC in G 359.1—0.5 are nearly one order of magnitude 
higher than the underlying bremsstrahlung (see Model A 
in figure 3), and kT^ determined from the width of the 
RRC structures is significantly away from kTciE deter- 
mined mainly from the bremsstrahlung continuum in the 
1-CIE plasma. Therefore, such a simple analysis cannot 
give physical quantities properly, such as total energy and 
metal abundances. 

We apply an over-ionized plasma model^ for the first 
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Fig. 3. Same spectrum as figure 2a, but fitted with Model A. 
The CIE plasma model, the RRC, and the Lya lines of the Si 
and S are shown as red, green, and blue lines, respectively. 
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Fig. 4. Same spectrum as figure 2a, but fitted with Model B. 

time to the whole SNR spectrum, in which electron and 
ionization temperatures are treated as independent free 
parameters (hereafter. Model B). The abundances of Mg, 
Si, and S are free parameters, while the other elements 
are fixed at the solar values. Model B gives the best-fit 
X^/d.o.f. of 223/189, which is marginally acceptable (null 
hypothesis probability = 0.05). The best-fit model and 
parameters of Model B are shown in figure 4 and table 3, 
respectively. 

3.2.4- Spatial structure of the spectra 

Although the size of G 359.1—0.5 is reasonably large 
compared to the angular resolution of the XIS, the poor 
photon statistics hmit the study of the spatial structure 
of the over-ionized plasma. We therefore simply divide 
the source region of G 359.1—0.5 into two at 6 = — 0.°5: 
the near side of the Galactic plane (here, near side) and 
the far side from the plane (far side), and examine the 



For this analysis, we used a code named "CIE" in the SPEX 



package. Although the name of "CIE" is misleading, this code 
can treat the electron (kTe) and ionization (kTz) temperatures 
as independent free parameters, and hence can describe plasma 
states in over-ionization {kTz > kTe). 
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Table 3. The best-fit parameters of Model B. 



Component 


Parameter 


Value 


Absorption 
Over-ionized plasma 


TVh (10^^ H cm-^) 


"^'^ -0.2 




kTe (keV) 


0.29 ±0.02 




kT^ (keV) 


r| 77 +0.09 
1 1 -0.08 




Zug (solar) 


o 4 +1.2 
'^•^ -0.9 




Zsi (solar) 


12 tl 




Zs (solar) 


17 tl 




VEM* 


ij 9 +2.3 
'•^ -1.8 



* Volume emission measure. The unit is j nen^dV / AttD^ 
(^10^^ cm~^) , where rie, rip, V ^ and D are the electron and proton 
densities (cm""^), the emitting volume (cm"^), and the distance to 
the source (cm), respectively. 



spectral difference. We extract the spectra from the near 
side and the far side of the source region separately, and fit 
each spectrum with the same Model B of the fixed column 
density (A/'h) and the abundances of Mg, Si, and S (^Mg, 
Zsi, and Zs) at the best-fit value of Model B (table 3). 

The best-fit parameters of the far side are kT^ of 0.26 
±0.03 keV, kT^ of 0.80 t^^W keV, and volume emission 
measure (VEM) of 3.0 ±0.3 xlO^^ cm-^ On the other 
hand, those of the near side are kT^ of 0.31 ±0.02 keV, 
kT^ of 0.75 tl^Q keV, and VEM of 4.7 ±0.4 xlO^^ cm-^ 

Although within the statistical errors, we see hints that 
the electron temperature in the far side is smaller than 
that in the near side and vice- versa for the ionization tem- 
perature. For the quantitative comparison, we define the 
ionization parameter (IP) as the ratio of ionization tem- 
perature to the electron temperature (/cTz/ZcTg), then we 
see that the IP in the far side of 3.0 is slightly higher 
than that in the near side of 2.4 ±0.1. From the best-fit 
VEM, we confirm that the density in the far side is smaller 
than that in the near side. 

Since the near side may be interacting with the Snake 
(figure 5), the non- thermal radio filament, the IP of this 
region might be influenced by high energy electrons, which 
may emit non-thermal X-rays by synchrotron radiation. 
Although we find no hint of excess X-rays in the hard 
X-ray 4.0-8.0 keV band as shown in figure 5, we try to 
estimate the upper-limit of the hard X-ray flux. 

We make the X-ray spectrum from the pentagonal re- 
gion (magenta), where the background is taken from the 
adjacent region shown by the dotted magenta lines in fig- 
ure 5. Since the statistics is limited, we fit the spectrum 
with a power-law model of fixed index 2.0 affected by pho- 
toelectric absorption of the fixed A^h = 2.0 x 10^^ H cm~^, 
the best-fit parameter of Model B. Then we constrain the 
flux upper-limit to be 2xl0~^^ ergs cm~^ s~^ at the 90% 
confidence level. 

Using this upper-limit of the X-ray flux, we estimate 
the energy of the non-thermal electrons. Since the spec- 
tral index and the flux derived from 1446 and 4790 MHz 
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Fig. 5. XIS close-up image in the 4.0-8.0 keV band. To 
avoid contamination by the emission from the calibration 
source, we excluded the 5.7-6.7 keV band. The radio map 
by the VLA 20 cm observation is superimposed in the con- 
tours (Yusef-Zadeh et al. 2004). 

vary along the length of the Snake, we adopt the spec- 
tral index to be 0, based on figure 12 of Gray et al. 
(1995) within latitude from —15' to —21', which corre- 
sponds to the pentagonal region. Assuming that the size 
of the Snake in the radio observations within the pentag- 
onal region is 6' x 0.'2, we estimate the flux density to 
be 70 mJy at 1446 MHz. We adopt the magnetic field 
strength of 88 /iG, which is the minimum magnetic field 
of the Snake estimated by the radio observations (Gray 
et al. 1995). Then the best-fit srcut model, which rep- 
resents the synchrotron spectrum from an exponentially 
cut-off power-law distribution of electrons (Reynolds 1998; 
Reynolds & Keohane 1999), gives the upper-limit of the 
cut-off electron energy to be 1 GeV. 

4. Discussion 

We find that the peculiar X-ray spectrum of 
G 359.1—0.5 is nicely reproduced by the over-ionized 
plasma model with kT^ of 0.29 and kT^ of 0.77 keV unlike 
the previous X-ray study (Bamba et al. 2000), in which the 
X-ray spectra required the two-component CIE plasma. 
Due to the limited energy resolution and photon statis- 
tics of ASCA, the RRC of Si was misunderstood as the 
Lya line of S, and thus the non-existent high-temperature 
(4.4 keV) S-rich plasma would be claimed. In the follow- 
ing sections, we discuss possible implications based on this 
over-ionized plasma model. 

4-1' Physical parameters in the over-ionized model 

We note that the Si and S abundances obtained from 
the Model A are different from those of Model B. This 
would demonstrate the limitations of simple analysis of 
Model A in high IP cases. 

We must also pay attention to the result of Model B that 
continuum X-rays due to the RRC from heavy elements 
are significant compared to the bremsstrahlung. Since 
the RRC emission in 1.0-2.0 keV is mainly due to the 
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abundant elements, O and Ne, the bremsstrahlung flux 
from the H and He atoms is greatly affected by the O and 
Ne abundances. We therefore cannot properly estimate 
the metal abundances relative to H. 

In order to estimate the abundance errors due to pu- 
tative large RRC of O and Ne, we fit the spectra with 
Model B, but the abundances of O and Ne are fixed as 
4 solar, near the same value of Mg. The apparent abun- 
dances of Mg, Si, and S are 9.9, 34, and 48 solar (x^/d.o.f. 
of 227/189), which are largely different from table 3 (for 
1-solar O and Ne). On the other hand if we estimate the 
abundances of Mg, Si, and S relative to O or Ne, then the 
relative abundances are modified to 2.5, 8.6, and 12 solar, 
for Mg, Si, and S, respectively. These values are similar 
to those of Model B (table 3). 

With the same reason, the VEM calculated from the 
bremsstrahlung of H and He is 2.4 x 10^^ cm~^ in 4-solar O 
and Ne abundances, which is about a third of the Model B. 
On the other hand, TVh, /cTe, and kT^ are estimated to be 
2.0 xlO^^ H cm"^ 0.29, and 0.79 keV, respectively, thus 
the uncertainty of these parameters are small. 

The best-fit column density of 2.0x10^^ H cm~^ (table 
3) is equal to or slightly smaller than that of other X-ray 
sources near to the GC with the same Galactic coordi- 
nates of this SNR (Sakano et al. 1999). Hence, this SNR 
would be located in or front of the GC. This is consistent 
with the previous radio observations. Hereafter, we as- 
sume that the SNR is located in the GC region with the 
distance of 8.5 kpc. The shape of the source region is ap- 
proximated with an ellipse of 12.0x19.6 pc^ (see figure 1). 
Assuming that the total kinetic energy of ions (protons) 
in the plasma is equal to /cTg, the sound velocity is esti- 
mated to be 3x10^ cm s~^. Dividing the major axis by 
the sound velocity, we estimate the dynamical time scale 
to be 7x10"^ yr. 

4-2. Origin of the over-ionized plasma 

The most remarkable discovery is the extremely 
over-ionized plasma from the mixed-morphology SNR 
G 359.1—0.5. The ionization temperature of kT^ = 
0.77 keV is similar to, but electron temperature of kT^ = 
0.29 keV is far lower than those of the middle-aged SNR. 
Although the number of samples for over-ionized SNRs 
are very small, G 359.1—0.5 is unusual compared to the 
other mixed morphology SNRs, IC443 and W49B. The 
possible scenarios for the origin of over-ionized plasmas in 
IC443 and W49B have been discussed (Kawasaki et al. 
2002; Yamaguchi et al. 2009; Miceh et al. 2006; Miceh 
et al. 2010). Unfortunately, the observational results of 
G 359.1—0.5 are very limited due to the poor statistics for 
the quantitative discussion on the origin, we hence briefly 
comment and address some plausible scenarios. 

Thermal conduction: Kawasaki et al. (2002) proposed 
that thermal conduction from a central hot plasma to a 
surrounding cool plasma caused the strong cooling of the 
electrons in IC443. G 359.1—0.5 was reported to be sur- 
rounded by the ^^CO molecular ring and the Hi emis- 
sion ring (Uchida et al. 1992a; Uchida et al. 1992b). 
Furthermore, OH maser (1720 MHz) sources, tracers of 



shock activity, are located at the projected position of the 
SNR sheh (Yusef-Zadeh et al. 1995; Yusef-Zadeh et al. 
1996). Since the emissions of the surrounding ^^CO, Hi, 
and OH masers are more intense at the shell of the near 
side than that of the far side, this part of G 359.1—0.5 
would interact strongly with the molecular cloud. Then 
the thermal conduction from the plasma to the surround- 
ing molecular cloud might cause the cooling of the elec- 
trons more effectively and produce a more over-ionized 
plasma in the near side. Thus, the IP in the near side 
may be higher than that in the far side. 

Adiabatic cooling: Yamaguchi et al. (2009) proposed 
that the drastic adiabatic expansion of the shocked gas 
causes the rapid cooling of electrons (Itoh & Masai 1989). 
The supernova explodes in a dense circumstellar medium, 
and then the gas is shock-heated and significantly ionized 
at the initial phase of the SNR evolution. The subsequent 
outbreak of the blast wave to a low-density interstellar 
medium (ISM) causes the drastic adiabatic expansion of 
the shocked gas and the resultant rapid cooling of the 
electrons. This scenario predicts that the IP is higher in 
the region with lower electron density (far side), since the 
adiabatic expansion in the low-density ISM produces an 
over-ionized plasma. 

Suprathermal electrons: Over-ionized plasmas were 
found in big solar flares with hard X-ray tails (e.g. Tanaka 
1986). Kato & Masai (1992) proposed that a few per- 
cent of the suprathermal electrons, which produce the 
hard X-rays by bremsstrahlung, caused effective ioniza- 
tion. Suprathermal electrons may be supplied by the ra- 
dio filament Snake. Since the most effective ionizations 
are made near a few 10 keV (by electrons) and several 
10 MeV (by protons), our constraint of the upper-limit of 
electron energy of 1 GeV does not rule-out this possibility. 

We find a hint that the IP is larger in the far side than 
in the near side. This fact may favor the second scenario, 
the adiabatic cooling of electrons. For quantitative calcu- 
lation on the origin of this extremely over-ionized plasma, 
detailed observations of the spatial structure is encour- 
aged. 
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